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ENVIRONMENTAL  OCEAN-BOTTOM  BACKSCATTERING  STRENGTH 
BASED  ON  THE  MOURAD-JACKSON  MODEL 


INTRODUCTION 

Scattering  in  the  ocean  subbottom  volume  and  scattering  by  the  interface  between  the  water  column 
and  the  subbottom  sediment  region  are  both  major  contributors  to  acoustic  reverberation  in  the  ocean  and 
thus  significantly  influence  the  performance  of  active  acoustic  systems.  This  is  especially  true  in  the  shallow 
water  that  describes  most  of  the  environments  present  in  littoral  (near-coastal)  regions.  In  these  regions,  it  is 
very  difficult  to  propagate  acoustic  energy  without  interacting  with  the  ocean  bottom. 

Most  acoustic  reverberation  and  active-acoustic-system  performance  prediction  models  (see,  for  ex¬ 
ample  Refs.  1,  2,  and  3)  used  by  the  Navy,  even  in  scientific  contexts,  incorporate  scattering  formulations 
that  do  not  include  physical  descriptors  of  the  ocean  subbottom  and  bottom-ocean  interface.  They  are  often 
either  based  on  ad  hoc  formulations  such  as  Lambert’s  Law  or  the  so-called  Urick-MacKenzie  [4,5]  ap¬ 
proach.  Although  such  formulations  have  been  successfully  used  in  experimental  planning  and  modeling 
experimental  results,  they  do  not  directly  incorporate  the  salient  environmental  features  that  apparently 
influence  scattering. 

Pierre  Mourad  and  Darrell  Jackson  of  the  Applied  Physics  Laboratory  at  The  University  of  Washington 
(APLAJW)  have  developed  a  formulation  [6]  for  bottom  backscattering  that  incorporates  both  scattering 
from  the  water-sediment  interface  and  scattering  from  the  subbottom  sediment  volume.  The  environmental 
contribution  to  backscattering  in  the  Mourad-Jackson  model  is  contained  in  eight  parameters.  These  param¬ 
eters  are  summarized  in  Table  1 . 

The  meaning  of  most  of  the  Mourad  and  Jackson’s  parameters  should  be  clear  from  their  definitions. 
The  loss  parameter  governs  energy  lost  while  propagating  through  the  sediment  volume.  The  volume  pa¬ 
rameter  governs  the  amount  of  scattering  in  the  sediment.  The  spectral  strength  and  spectral  exponent  are 
measures  of  the  roughness  of  the  water-sediment  interface.  All  parameters  except  water  sound  speed  (m/s), 
sound-speed  gradient  (1/s),  and  spectral  strength  (m'')  are  dimensionless. 

The  complete  Mourad-Jackson  formulation,  in  the  form  that  we  considered  it,  is  described  in  detail  in 
Ref.  6  and  its  references.  Besides  for  these  environmental  parameters,  the  backscatter  also  depends  on  the 
single  acoustic  frequency /of  the  source  and  the  incident  or  grazing  angle  0  (measured  from  the  horizontal) 
of  the  energy  interacting  with  the  boundary. 

The  eight  environmental  parameters,  although  mathematically  independent,  are  not  unrelated.  For  ex¬ 
ample,  the  speed  of  sound  in  the  water  at  the  water-sediment  interface  appears  in  the  sound-speed  ratio  v, 
besides  for  being  c^.  Also,  relationships  that  express  one  parameter  in  terms  of  a  perhaps  more  fundamental 
quantity  and  other  parameters  are  sometimes  postulated.  One  example  of  this  is  the  following  equation. 
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Table  1-  Environmental  Parameters 


Symbol 

Definition 

Short  Name 

S 

Sound  speed  in  water  at  water/sediment 
interface  (m/s) 

Water  sound  speed 

P 

Ratio  of  sediment  mass  density  to  water  mass 
density 

Density  ratio 

V 

Ratio  of  sediment  sound  speed  (at 
water/sediment  interface)  to  water  sound 
speed 

Sound-speed  ratio 

g 

Slope  of  sediment  sound-speed  profile  at 
water/sediment  interface  (s'') 

Sound-speed  gradient 

8 

Ratio  of  imaginary  wave  number  to  real  wave 
number  for  the  sediment 

Loss  parameter 

Ratio  of  sediment  volume  scattering  cross 
section  to  sediment  attenuation  coefficient 

Volume  parameter 

Y 

Exponent  of  bottom  relief  spectrum 

Spectral  exponent 

w. 

Strength  of  bottom  relief  spectram  (m'*)  at 
wave  number  1  rad/m 

Spectral  strength 

which  expresses  the  loss  parameter  in  terms  of  other  parameters,  and  a^,  which  is  the  attenuation  coefficient 
in  dB/m  (assuming  that  is  in  m/s)  [7], 


(XnvCf)  In  10 

a) 

The  remainder  of  such  relationships  are  defined  and  described  by  Mourad  and  Jackson.  A  positive  value  of 
the  sound-speed  gradient  g  indicates  sound  speed  increases  with  depth  in  the  sediment. 

The  long-term  goal  of  this  research  is  to  develop  an  understanding  of  the  spatial  and  temporal  charac¬ 
teristics  of  the  shallow-water  monostatic  and  bistatic  reverberation  fields  by  investigating  the  dominant 
environmental  features  as  a  function  of  the  acoustic  frequency  and  the  source-receiver  configuration.  More 
specifically,  the  objectives  of  this  report  are  to  study  the  dominant  environmental  parameters  governing 
ocean  bottom  backscattering  using  the  Mourad- Jackson  formulation  and  decide  on  the  appropriateness  of 
the  formulation  as  a  component  of  an  active-acoustic-system  performance  prediction  model. 

A  scattering  strength  formulation  using  this  model  for  backscattering  strength  is  being  considered  as 
an  alternative  for  Lambert’s  Law  in  the  BiRASP  [2]  and  BiKR  [8]  active-acoustic-system  performance 
prediction  models.  The  Mourad- Jackson  formulation  results  in  direct  backscatter  strength,  which  is  a  mea- 
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sure  of  the  energy  backscattered  at  the  same  angle  0  as  the  incoming  energy.  To  be  useful  for  BiRASP  or 
BiKR,  assumptions  must  be  made  as  to  how  the  energy  varies  in  angles  other  than  the  backscatter  direction. 
Such  assumptions  are  not  discussed  here. 

The  direct  backscatter  strength  predicted  by  the  Mourad-Jackson  formulation  is  the  sum  of  two  terms, 
one  representing  the  energy  scattered  by  the  irregular  bathymetry,  and  the  other  the  energy  scattered  by 
random  inhomogeneities  in  the  subbottom  volume.  The  scattering,  in  turn,  depends  on  the  parameters  listed 
in  Table  1. 

Originally,  we  attempted  to  analytically  study  how  the  backscattered  field  is  influenced  by  each  of 
Mourad  and  Jackson’s  parameters.  We  examined  the  water-sediment  interface  backscatter  and  the  subbottom 
volume  backscatter,  both  individually  and  together.  Although  a  few  general  statements  were  learned,  it 
eventually  became  obvious  that  an  analytic  approach  would  not  lead  to  general  understanding  of  the  scatter¬ 
ing  mechanisms.  This  is  because  Mourad  and  Jackson  had  already  introduced  computational  simplifications 
that  would  not  allow  at  this  point  further  analytic  simplifications.  Their  formulation  resulted  in  a  complex 
combination  of  transcendental  functions  and  summations  that,  although  quite  computationally  feasible,  does 
not  lend  itself  to  any  analytic  simplification.  Thus,  further  work  using  the  analytic  approach  was  not  contin¬ 
ued  and  is  not  included  here. 

The  most  obvious  fact  learned  by  examining  the  equations  in  Mourad  and  Jackson’s  formulation  is 
that  the  volume  parameter  only  influences  the  bottom  volume  backscatter.  It  does  not  contribute  to  the 
interface  backscatter  calculation. 

NUMERICAL  RESULTS 

When  an  analytic  approach  proved  to  be  impractical,  we  decided  to  use  a  numerical  approach.  We 
chose  to  calculate  both  interface  scatter  and  bottom  volume  contributions  to  the  direct  backscatter  strength 
numerically.  We  did  this  by  using  ranges  for  each  parameter  that  span  the  values  expected  to  be  found  in  the 
ocean,  with  emphasis  on  littoral  regions.  Such  parameter  ranges  were  obtained  by  us  both  directly  from 
cases  used  by  Mourad  and  Jackson  in  their  published  work  and  from  a  substantial  study  [9]  of  parameter 
values  by  Anthony  Lyons,  then  of  Texas  A  &  M  University,  using  research  of  Hamilton  [10],  Hamilton  and 
Bachman  [11],  and  others  (including  Mourad  and  Jackson’s  own  work). 

The  numerical  approach  proved  to  be  more  productive  than  the  analytic  approach.  However,  it  is  still 
not  practical  with  existing  computers  and  workstations  to  calculate  the  backscattering  strength  for  all  values 
of  the  10-parameter  set  of  values  formed  from  combining  the  8  environmental  parameters  with  the  acoustic 
frequency /and  the  incident  angle  0.  However,  preliminary  studies  of  variations  in  the  backscatter  function 
as  the  parameters  varied  over  their  ranges  that  Lyons  produced  showed  that  we  could  safely  leave  some  of 
the  parameters  fixed. 

We  thus  chose  a  constant  value  of  =  1530  m/s.  Although  the  backscattering  strength  depends  on  the 
water  sound  speed  at  the  water-sediment  interface,  the  variation  is  quite  smooth.  Also,  changes  in  the  sound 
speed  can  produce  similar  results  to  corresponding  changes  in  the  acoustic  frequency,  and  such  changes  are 
analyzed  here.  Thus,  major  variations  in  backscatter  strength  are  not  omitted  by  taking  a  constant  value  for 
the  parameter  c^. 

As  observed,  the  volume  parameter  does  not  enter  into  the  interface  backscattering  calculation,  and 
it  enters  into  the  subbottom  volume  backscattering  calculation  only  as  a  scaling  parameter.  Thus,  with  no 
loss  of  generality  in  the  calculations,  we  pick  =  0.0003  as  a  typical  value. 
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Lastly,  little  is  known  from  direct  at-sea  measurements  about  the  spectral  exponent  y.  Mourad  and 
Jackson  have,  in  general,  assumed  a  value  of  Y=  3.25  for  this  parameter  because  of  this  lack  of  experimental 
knowledge.  We,  therefore,  also  assumed  this  value  for  all  cases  reported  here. 

There  remain  five  environmental  parameters  in  addition  to  /and  0.  We  allowed  the  incident  (grazing) 
angle  0  to  vary  from  0°  to  60°  in  steps  of  2°  (31  values).  We  used  acoustic  frequencies  of  50,  100,  200, 400, 
800,  and  1600  Hz,  which  span  the  low-  to  mid-frequency  range  of  interest  in  many  past  active  acoustic 
experiments  performed  by  NRL  acousticians. 

We  desired  to  study  the  behavior  of  the  backscattering  strength  produced  by  the  Mourad-Jackson  model 
mainly  in  the  form  of  color  or  gray-scale  plots  with  0  as  the  abscissa  and  one  of  the  other  five  parameters  (p, 
V,  5 ,  g,  and  w^)  as  the  ordinate.  Therefore,  we  required  fairly  detailed  sampling  (similar  to  the  sampling  used 
for  angle)  for’the  parameter  used  as  the  ordinate  and  coarser  sampling  for  the  remaining  parameters,  which 
were  held  constant  in  a  given  plot.  We  thus  chose  parameter  values  as  listed  in  Table  2. 


Table  2  -  Parameter  Values 


Parameter 

Values  when  used  as  y-axis 
start  (increment)  end 

Parameter  values  used  for  individual 

runs 

P 

1.0  (0.05)  2.5 

1.00,  1.25,  1.50,  2.00,  2.50 

V 

0.9  (0.01)  1.3 

0.90,  0.95,  1.00,  1.10,  1.25,  1.50 

6 

0.0001  (0.001)  0.0301 

0.0001,  0.001,  0.01,  0.02,  0.03 

8 

-0.1  (-0.5)  -18.1 

-0.1,  -0.5,  -1.0,  -10.0,  -18.0 

0.00001  (0.005)  0.18001 

0.00001,  0.0001,  0.001,  0.01,  0.1 

The  parameter  used  as  the  ordinate  is  sampled  by  a  set  of  3 1  to  41  values,  as  indicated  in  the  second 
column  of  Table  2.  The  three  numbers  in  this  column  are  respectively,  the  lowest  value  of  the  parameter 
considered,  the  increment  in  value  of  the  parameter,  and  the  largest  value  of  the  parameter  used.  The  data¬ 
base  was  completed  by  using  the  set  of  five  or  six  values  of  the  remaining  four  parameters,  as  indicated  in 
the  third  column.  For  example,  we  generated  a  set  of 4,500  files  to  generate  plots  where  density  ratio  p  is  the 
ordinate,  corresponding  to  the  4,500  combinations  of  the  six  frequencies,  six  values  of  the  sound-speed 
ratio,  five  values  of  the  loss  parameter,  five  values  of  the  sound-speed  gradient,  and  five  values  of  the 
spectral  strength.  Similar  sets  were  generated  with  each  of  the  other  four  environmental  parameters  as  the 
ordinate.  There  are  4,500  files  for  each  choice  of  parameter  for  the  ordinate,  except  for  only  3,750  when  the 
sound-speed  ratio  v  is  the  ordinate,  for  a  total  of  21,750  files. 

It  is  possible  that  some  combinations  of  parameter  values  are  physically  unobtainable.  That  is,  no 
region  of  the  ocean  floor  is  constructed  of  a  material  that  results  in  certain  combinations  of  parameter 
values.  Thus,  by  using  the  set  of  parameter  values  described  in  Table  2,  we  expect  to  have  more  than  covered 
the  region  of  parameter  values  that  exist  in  the  world’s  littoral  regions. 
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Note  that  separate  runs  of  the  Mourad-Jackson  model,  as  originally  implemented,  correspond  to  each 
grazing  angle,  so  this  database  represents  more  than  2.3  x  10’  model  mns.  Because  of  the  large  amount  of 
computer  time  needed  to  generate  the  entire  database,  about  a  half-dozen  Silicon  Graphics  processors  were 
each  assigned  part  of  the  task.  It  still  took  several  weeks  of  wall  clock  time  to  produce  the  database. 

The  computational  time  for  the  21,750  files  varied  by  a  factor  of  more  than  three  orders  of  magnitude 
between  slowest  and  fastest,  depending  mainly  on  the  time  needed  for  numerical  integration  of  the  subbottom 
backscatter  to  converge.  The  water-sediment  interface  backscattering  contribution  required  virtually  no  time 
to  calculate  compared  to  the  subbottom  backscatter.  If  computational  times  for  the  individual  runs  were 
compared,  the  times  would  vary  by  even  much  more  than  three  orders  of  magnitude. 

Almost  universally,  higher  frequency  calculations  took  more  central  processor  time  than  lower  fre¬ 
quency  calculations.  Also,  calculation  time  depended  greatly  on  the  sound-speed  gradient  g  -  the  larger 
gradient  calculations  being  much  faster  than  the  calculations  with  smaller  gradients.  This  is  an  indication  of 
deeper  penetration  into  the  sediment  when  the  gradient  is  smaller.  Integrations  used  in  computing  the 
subbottom  backscatter  take  much  longer  to  converge  with  deeper  penetration. 

It  is  impossible  to  illustrate  this  report  with  plots  corresponding  to  each  of  the  21,750  files.  Thus,  only 
selective  files  are  shown  here.  However,  all  of  the  files  not  contained  in  this  report  have  been  examined  on 
Silicon  Graphics  workstation  monitors. 

Water-Sediment  Interface  Backscatter 

Figure  1  is  a  plot  showing  variation  of  the  surface  component  of  the  backscatter  strength  as  a  function 
of  the  sound-speed  ratio  v  and  grazing  angle  0.  The  fixed  parameters  are  the  source  frequency/=  100  Hz,  the 
sound-speed  gradient  g  =  0.1  s  ',  the  density  ratio  p  =  1.0,  the  loss  parameter  5  =  0.0001,  and  the  spectral 
strength  =  0. 1  m'*.  The  spectral  strength  value  is  quite  high,  although  not  unrealistic.  The  color  scale  as 
shown  varies  from -20  dB  to  -100  dB.  We  note  a  generally  smooth,  slow  variation  in  the  backscatter  strength 
except  for  a  null  of  at  least  60  dB  at  v  =  1.0.  This  value  of  the  sound-speed  ratio,  combined  with  the  fixed 
density  ratio  p  =  1 .0,  essentially  corresponds  to  there  being  no  interface  (it  is  acoustically  “invisible”).  Thus, 
one  would  expect  to  see  minimal  backscattering  strength  from  the  surface,  even  though  it  is  rough.  Also,  as 
might  be  expected,  there  is  no  surface  component  to  the  backscatter  at  0°  grazing  angle. 

Besides  for  the  variations  of  the  types  discussed  in  the  previous  paragraph,  the  surface  component  of 
the  backscatter  strength  exhibits  little  structure  over  the  ranges  of  values  considered.  Thus,  further  plots  of 
this  type  are  not  shown.  However,  we  note  that  the  surface  component,  although  not  having  much  structure, 
can,  for  appropriate  choices  of  parameter  values,  be  much  greater  in  value  than  the  subbottom  component  of 
the  backscatter  strength.  In  addition,  adjusting  the  value  of  the  volume  parameter  will  change  the  relation¬ 
ship  between  the  interface  and  the  subbottom  contributions  to  bottom  backscatter  as  predicted  by  the  Mourad- 
Jackson  formulation  because,  as  already  stated,  this  parameter  enters  into  the  backscatter  calculation  strictly 
as  a  scaling  factor. 

Subbottom  Backscatter  -  Typical  Parameter  Values 

In  the  next  five  plots  (Figs.  2-6),  we  investigate  variations  in  the  subbottom  component  of  the  backscat¬ 
ter  strength.  We  do  this  by  using  typical,  midrange  values  for  each  of  the  environmental  parameters  that  is 
being  held  fixed,  and  allowing  each  parameter  to  be  used,  in  turn,  as  the  ordinate.  The  environmental  param¬ 
eters  are  taken  as  sound-speed  gradient  g  =  1.0  s  ',  sound-speed  ratio  v  =  1.1,  density  ratio  p  =  1.5,  loss 
parameter  5  =  0.01,  and  spectral  strength  =  0.001  m**.  We  also  choose  a  source  frequency  in  the  middle  of 
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Fig.  1  -  Interface  backscatter  strength  as  function  of  sound>speed  ratio  v  and  grazing  angle  0.  Frequency /=  100  Hz. 
Environmental  parameters:  g  =  0.1  s'*,  p  =  1, 5  =  0.0001,  and  =  0.1 


dB  Level 


Fig.  2  -  Sediment  backscatter  strength  as  function  of  density  ratio  p  and  grazing  angle  0.  Frequency/^  400  Hz.  Environmental 
parameters:  g  =  1.0  s'^  v  =  1.1,  p  =  0.01,  and  =  0.001  m"*.  Mark  on  vertical  axis  indicates  value  common  to  Figs.  2-6. 
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the  range  of  values  treated,  /=  400  Hz.  The  small  mark  on  the  y-axis  corresponds  to  the  value  chosen  for  that 
parameter  in  the  remaining  four  plots. 

First,  in  Fig.  2,  we  allow  the  density  ratio  to  vary  and  the  other  parameters  to  be  as  in  the  previous 
paragraph.  The  subbottom  backscattering  strength  exhibits  smooth  variation  as  a  function  of  density  ratio 
and  grazing  angles  for  the  remainder  of  the  “typical”  parameter  values.  There  is  very  little  backscatter  at 
zero  degrees  grazing  because  of  lack  of  bottom  penetration,  and  the  backscatter  increases  monotonically 
with  increasing  grazing  angle.  At  low  grazing  angles,  the  backscatter  increases  with  increasing  density  ratio, 
while  at  higher  grazing  angle  it  decreases;  however  all  of  these  variations  are  quite  smooth.  The  smoothness 
of  the  variations  indicates  that  the  behavior  could  be  approximated  by  interpolating  between  a  relatively 
few,  well  chosen  data  points.  Thus,  the  formulation  appears,  for  this  set  of  parameter  values,  to  be  usable  as 
a  component  of  an  active-acoustic-system  performance  prediction  model. 

In  Fig.  3,  we  allow  the  strength  of  the  spectrum  of  the  bottom  relief  (spectral  strength)  to  vary.  Note  that 
the  typical  value  of  this  parameter,  used  in  Figs.  2  and  4-6,  is  near  the  x-axis.  This  is  because  of  the  great 
variation  in  magnitude  of  the  observed  values  of  the  exponent  in  the  bottom  relief  spectrum.  As  in  the 
previous  figure,  the  variations  are  quite  smooth,  indicating  that  the  Mourad  Jackson  model  is  suitable  for 
systems  performance  prediction. 


In  Fig.  4,  we  allow  the  sound-speed  gradient  in  the  sediment  at  the  water-sediment  interface  to  vary. 
There  is  somewhat  more  structure  here  than  in  either  of  the  two  previous  plots,  especially  at  grazing  angles 
greater  than  about  30  .  However,  even  here,  the  variations  appear  to  be  relatively  minor,  in  comparison  to 
the  magnitude  of  the  subbottom  component  of  the  backscattering  strength  itself,  thus  not  precluding  a  simple 
fit  to  the  data.  Also,  some  of  the  apparent  structure  in  the  upper  right  of  this  plot  is  due  to  a  slight  under¬ 
sampling  because  there  are  only  1147  data  points  (31  x  37).  A  more  severe  case  of  undersampling  is  illus¬ 
trated  in  the  next  section  of  this  report.  This  region  actually  contains  a  low-magnitude  periodic  variation, 
which  would  become  more  apparent  at  lower  frequency. 

In  Fig.  5,  we  vary  the  ratio  of  the  sediment  sound  speed  to  the  water  sound  speed  at  the  interface.  For 
the  range  of  values  of  v  considered,  the  subbottom  component  of  the  backscattering  strength  has  a  maximum 
value  for  grazing  angles  less  than  about  40°.  This  corresponds  approximately  to  v  =  1.0  for  0  =  0  and 
increases  with  increasing  grazing  angle.  The  backscattering  strength  increases  with  increasing  grazing  angle. 
Note,  however,  that  near  the  point  of  the  maximum  for  vertical  cuts,  the  rate  of  increase  with  increasing 
grazing  angles  decreases. 


Lastly,  in  Fig.  6,  we  allow  the  ratio  of  imaginary  wave  number  to  real  wave  number  in  the  sediment  (the 
loss  parameter)  to  vary.  As  in  the  previous  four  cases,  the  subbottom  component  of  the  backscattering  strength 
varies  fairly  smoothly.  At  a  grazing  angle  of  approximately  20°,  the  backscatter  strength  changes  from 
increasing  with  increasing  loss  parameter  to  decreasing  with  increasing  loss  parameter.  There  is  a  null  for 
5  =  0. 

In  summary.  Figs.  2-6  illustrate  that  for  typical  values  of  the  environmental  parameters,  the  subbottom 
component  of  the  backscattering  strength  does  not  have  any  extreme  variations.  Therefore,  the  use  of  the 
Mourad- Jackson  formulation  as  a  component  of  an  active-acoustic-system  performance  prediction  model  is 
not  ruled  out  from  considerations  of  computational  efficiency.  Most  of  the  more  than  20,000  situations 
studied  look  similar  to  the  corresponding  plot  in  Figs.  2-6,  and  so  this  same  conclusion  should  hold  in  those 
cases  also. 
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3  -  Sediment  backscatter  strength  as  function  of  spectral  strength  and  grazing  angle  0  .  Frequency /=  400  Hz.  Environ¬ 
mental  parameters:  g  =  1.0  s’’,  v  =  1.1,  p  =  1.5,  and  5  =  0.01.  Mark  on  vertical  axis  indicates  value  common  to  Figs.  2-6. 
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Fig.  4  -  Sediment  backscatter  strength  as  function  of  sound-speed  gradient  g  and  grazing  angle  0  .  Frequency/  =  400  Hz.  En¬ 
vironmental  parameters:  v  =  1.1,  p  =  1.5, 8  =  0.01,  and  =  0.001  m^  Mark  on  vertical  axis  indicates  value  common  to  Figs.  2-6. 
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Figures  7(a)  and  7(b)  are  analogous  to  Fig.  2,  with  the  density  ratio  being  allowed  to  vary.  However,  in 
Fig.  7(a),  the  frequency  is/=  50  Hz,  and  in  Fig.  7(b),  the  frequency  is/=  1600  Hz.  The  general  level  of  the 
plots  is  similar  to  the  original  400  Hz  case,  although,  of  course,  the  structure  is  not  identical. 

Similar  plots  to  Figs.  3-6  at  50  Hz  and  1600  Hz  would  also  have  similar  levels,  except  in  some  cases 
structural  features  would  be  more  pronounced.  For  example,  in  the  50  Hz  case  analogous  to  Fig.  4  (sound- 
speed  gradient  varying),  there  is  a  clearly  periodic  stracture  in  which  both  grazing  angle  and  sound-speed 
gradient  are  high.  The  entire  structure  is  missing  at  1600  Hz. 

We  now  turn  our  attention  to  some  situations  where  more  structure  is  apparent  in  the  subbottom  com¬ 
ponent  of  the  backscattering  strength  than  in  Figs.  2-6.  In  Figs.  8(a)  and  8(b)  (as  in  Fig.  3),  the  spectral 
strength  is  allowed  to  vary.  The  environmental  parameters  are  the  sound-speed  gradient  g  =  10.0  s'',  the 
sound-speed  ratio  v  =  1.25,  the  density  ratio  p  =  1.5,  and  the  loss  parameter  5  =  0.01.  In  Fig.  8(a),  the 
acoustic  frequency  is  50  Hz,  and  in  Fig.  8(b),  it  is  100  Hz.  The  values  of  the  density  ratio  and  the  loss 
parameter  are  as  in  the  previous  plots  while  the  sound-speed  gradient  is  greater  by  a  factor  of  10. 

There  is  significantly  more  stmcture  in  both  of  these  cases  than  in  the  early  figures.  We  note  a  quasi- 
periodic  pattern  from  the  upper  left  to  the  upper  right  in  each  of  the  plots,  with  twice  as  many  cycles  at  100 
Hz  as  at  50  Hz.  In  addition,  certain  grazing  angles,  which  correspond  to  the  angles  where  the  maxima  of  the 
sinusoid-like  structure  intersect  the  =  0  axis,  give  local  maxima  independent  of  w^.  Each  of  these  maxima 
is  more  than  10  dB  stronger  than  the  nearby  minimum. 

Although  the  complicated  structure  scene  of  the  subbottom  backscatter  illustrated  in  Figs.  8(a)  and  8(b) 
would  appear  to  preclude  the  use  of  a  simple  approximation  scheme,  this  is  probably  not  the  case.  The 
periodicity  of  the  stracture,  the  uncertainty  of  knowledge  of  environmental  parameter  values,  and  the  aver¬ 
aging  that  takes  place  in  active-acoustic-system  performance  calculations  is  expected  to  still  allow  for  ap¬ 
proximation  in  this  and  similar  cases  simply  by  smoothing  the  data,  assuming  that  the  stracture  is  known  and 
the  smoothing  is  performed  carefully. 

Subbottom  Backscatter  -  Atypical  Parameter  Values 

In  this  section,  we  consider  a  case  containing  more  extreme  values  in  the  parameter  space  although  still 
within  the  ranges  of  the  parameters  contained  in  Lyons’  study  [9]. 

In  Fig.  9(a),  the  environmental  parameters  chosen  are  the  sound-speed  ratio  v  =  1.0,  the  density  ratio 
p  =  2.0,  the  loss  parameter  5  =  0.0001,  and  the  spectral  strength  ^2=  0.0001  m‘*.  The  frequency  is  800  Hz. 
The  loss  parameter  is  very  small,  and  thus  there  is  little  attenuation  in  the  subbottom,  allowing  a  lot  of 
energy  to  return  to  the  ocean  volume.  Much  of  the  apparent  stracture  in  the  subbottom  backscatter  seen  in 
the  figure  is  due  to  undersampling  problems.  When  100  times  as  many  data  points  are  plotted  (ten  times  in 
each  direction).  Fig.  9(b)  results.  There  is  still  an  undersampling  problem  at  combinations  of  high  grazing 
angle  and  low  sound-speed  gradient,  but  the  pattern  is  now  clear.  Again,  even  though  there  is  visible  strac¬ 
ture,  the  use  of  properly  smoothed  data  will  not  preclude  usage  in  an  active-acoustic-system  performance 
prediction  model. 

CONCLUSIONS  AND  RECOMMENDATIONS 

We  have  examined  a  great  many  output  sets  of  the  Mourad-Jackson  scattering  formulation  that  cover 
the  span  of  values  indicated  by  Anthony  Lyons’  literature  search  [9].  A  small  subset  of  these,  consisting  of 
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Fig.  7  -  Sediment  backscatter  strength  as  function  of  density  ratio  p  and  grazing  angle  6 
Environmental  parameters:  g  =  1.0  s"',  v=  1.1,5  =  0.01,  and  w  =  0.001  m^. 
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both  typical  and  nontypical  cases,  illustrate  this  report.  There  is,  in  general,  very  little  structure  in  the  water- 
sediment  interface  contribution  to  the  backscatter.  All  variations  are  quite  benign. 

There  is  more  structure  in  the  subbottom  contribution  to  the  backscatter  when  some  of  the  environmen¬ 
tal  parameters  are  permitted  to  vary.  However,  for  nonpathological  choices  of  parameters,  the  variations 
appear  to  be  such  that  the  backscatter  field  can  be  approximated  by  interpolation  between  a  relatively  few 
well-chosen  data  points. 

It  is  thus  felt  that  the  Mourad-Jackson  formulation  represents  a  viable  candidate  for  incorporation  as 
part  of  a  scattering  kernel  in  active-acoustic-system  performance  prediction  models  such  as  BiKR  and  BiRASP. 
It  is  recommended  that  studies  be  initiated  in  how  to  do  this  efficiently.  Also,  studies  are  needed  to  determine 
a  method  for  efficiently  extending  such  approximations  into  a  general  scattering  kernel  as  opposed  to  the 
backscatter  contribution  discussed  in  this  report. 
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